Individuals with JAK1 variants are affected by syndromic features encompassing autoimmunity, atopy, colitis, and dermatitis.

Inborn errors of immunity lead to autoimmunity, inflammation, allergy, infection, and/or malignancy. Disease-causing JAK1 gain-of-function (GoF) mutations are considered exceedingly rare and have been identified in only four families. Here, we use forward and reverse genetics to identify 59 individuals harboring one of four heterozygous JAK1 variants. In vitro and ex vivo analysis of these variants revealed hyperactive baseline and cytokine-induced STAT phosphorylation and interferon-stimulated gene (ISG) levels compared with wild-type JAK1. A systematic review of electronic health records from the BioME Biobank revealed increased likelihood of clinical presentation with autoimmunity, atopy, colitis, and/or dermatitis in JAK1 variant-positive individuals. Finally, treatment of one affected patient with severe atopic dermatitis using the JAK1/JAK2-selective inhibitor, baricitinib, resulted in clinically significant improvement. These findings suggest that individually rare JAK1 GoF variants may underlie an emerging syndrome with more common presentations of autoimmune and inflammatory disease (JAACD syndrome). More broadly, individuals who present with such conditions may benefit from genetic testing for the presence of JAK1 GoF variants.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation.

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Unambiguous detection of SARS-CoV-2 subgenomic mRNAs with single-cell RNA sequencing.

Single-cell RNA sequencing (scRNA-Seq) studies have provided critical insight into the pathogenesis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19). scRNA-Seq library preparation methods and data processing workflows are generally designed for the detection and quantification of eukaryotic host mRNAs and not viral RNAs. Here, we compare different scRNA-Seq library preparation methods for their ability to quantify and detect SARS-CoV-2 RNAs with a focus on subgenomic mRNAs (sgmRNAs). We show that compared to 10X Genomics Chromium Next GEM Single Cell 3' (10X 3') libraries or 10X Genomics Chromium Next GEM Single Cell V(D)J (10X 5') libraries sequenced with standard read configurations, 10X 5' libraries sequenced with an extended length read 1 (R1) that covers both cell barcode and transcript sequence (termed "10X 5' with extended R1") increase the number of unambiguous reads spanning leader-sgmRNA junction sites. We further present a data processing workflow, single-cell coronavirus sequencing (scCoVseq), which quantifies reads unambiguously assigned to viral sgmRNAs or viral genomic RNA (gRNA). We find that combining 10X 5' with extended R1 library preparation/sequencing and scCoVseq data processing maximizes the number of viral UMIs per cell quantified by scRNA-Seq. Corresponding sgmRNA expression levels are highly correlated with expression in matched bulk RNA-Seq data sets quantified with established tools for SARS-CoV-2 analysis. Using this scRNA-Seq approach, we find that SARS-CoV-2 gene expression is highly correlated across individual infected cells, which suggests that the proportion of viral sgmRNAs remains generally consistent throughout infection. Taken together, these results and corresponding data processing workflow enable robust quantification of coronavirus sgmRNA expression at single-cell resolution, thereby supporting high-resolution studies of viral RNA processes in individual cells. IMPORTANCE Single-cell RNA sequencing (scRNA-Seq) has emerged as a valuable tool to study host-virus interactions, especially for coronavirus disease 2019 (COVID-19). Here we compare the performance of different scRNA-Seq library preparation methods and sequencing strategies to detect SARS-CoV-2 RNAs and develop a data processing workflow to quantify unambiguous sequence reads derived from SARS-CoV-2 genomic RNA and subgenomic mRNAs. After establishing a workflow that maximizes the detection of SARS-CoV-2 subgenomic mRNAs, we explore patterns of SARS-CoV-2 gene expression across cells with variable levels of total viral RNA, assess host gene expression differences between infected and bystander cells, and identify non-canonical and lowly abundant SARS-CoV-2 RNAs. The sequencing and data processing strategies developed here can enhance studies of coronavirus RNA biology at single-cell resolution and thereby contribute to our understanding of viral pathogenesis.

A comparative study of in vitro air-liquid interface culture models of the human airway epithelium evaluating cellular heterogeneity and gene expression at single cell resolution.

BACKGROUND: The airway epithelium is composed of diverse cell types with specialized functions that mediate homeostasis and protect against respiratory pathogens. Human airway epithelial (HAE) cultures at air-liquid interface are a physiologically relevant in vitro model of this heterogeneous tissue and have enabled numerous studies of airway disease. HAE cultures are classically derived from primary epithelial cells, the relatively limited passage capacity of which can limit experimental methods and study designs. BCi-NS1.1, a previously described and widely used basal cell line engineered to express hTERT, exhibits extended passage lifespan while retaining the capacity for differentiation to HAE. However, gene expression and innate immune function in BCi-NS1.1-derived versus primary-derived HAE cultures have not been fully characterized. METHODS: BCi-NS1.1-derived HAE cultures (n = 3 independent differentiations) and primary-derived HAE cultures (n = 3 distinct donors) were characterized by immunofluorescence and single cell RNA-Seq (scRNA-Seq). Innate immune functions were evaluated in response to interferon stimulation and to infection with viral and bacterial respiratory pathogens. RESULTS: We confirm at high resolution that BCi-NS1.1- and primary-derived HAE cultures are largely similar in morphology, cell type composition, and overall gene expression patterns. While we observed cell-type specific expression differences of several interferon stimulated genes in BCi-NS1.1-derived HAE cultures, we did not observe significant differences in susceptibility to infection with influenza A virus and Staphylococcus aureus. CONCLUSIONS: Taken together, our results further support BCi-NS1.1-derived HAE cultures as a valuable tool for the study of airway infectious disease.

A comparative study of in vitro air-liquid interface culture models of the human airway epithelium evaluating cellular heterogeneity and gene expression at single cell resolution.

The airway epithelium is composed of diverse cell types with specialized functions that mediate homeostasis and protect against respiratory pathogens. Human airway epithelial cultures at air-liquid interface (HAE) are a physiologically relevant in vitro model of this heterogeneous tissue, enabling numerous studies of airway disease 1â€"7 . HAE cultures are classically derived from primary epithelial cells, the relatively limited passage capacity of which can limit experimental methods and study designs. BCi-NS1.1, a previously described and widely used basal cell line engineered to express hTERT, exhibits extended passage lifespan while retaining capacity for differentiation to HAE 5 . However, gene expression and innate immune function in HAE derived from BCi-NS1.1 versus primary cells have not been fully characterized. Here, combining single cell RNA-Seq (scRNA-Seq), immunohistochemistry, and functional experimentation, we confirm at high resolution that BCi-NS1.1 and primary HAE cultures are largely similar in morphology, cell type composition, and overall transcriptional patterns. While we observed cell-type specific expression differences of several interferon stimulated genes in BCi-NS1.1 HAE cultures, we did not observe significant differences in susceptibility to infection with influenza A virus and Staphylococcus aureus . Taken together, our results further support BCi-NS1.1-derived HAE cultures as a valuable tool for the study of airway infectious disease.

Autoimmunity in Down's syndrome via cytokines, CD4 T cells and CD11c+ B cells.

Down's syndrome (DS) presents with a constellation of cardiac, neurocognitive and growth impairments. Individuals with DS are also prone to severe infections and autoimmunity including thyroiditis, type 1 diabetes, coeliac disease and alopecia areata1,2. Here, to investigate the mechanisms underlying autoimmune susceptibility, we mapped the soluble and cellular immune landscape of individuals with DS. We found a persistent elevation of up to 22 cytokines at steady state (at levels often exceeding those in patients with acute infection) and detected basal cellular activation: chronic IL-6 signalling in CD4 T cells and a high proportion of plasmablasts and CD11c+TbethighCD21low B cells (Tbet is also known as TBX21). This subset is known to be autoimmune-prone and displayed even greater autoreactive features in DS including receptors with fewer non-reference nucleotides and higher IGHV4-34 utilization. In vitro, incubation of naive B cells in the plasma of individuals with DS or with IL-6-activated T cells resulted in increased plasmablast differentiation compared with control plasma or unstimulated T cells, respectively. Finally, we detected 365 auto-antibodies in the plasma of individuals with DS, which targeted the gastrointestinal tract, the pancreas, the thyroid, the central nervous system, and the immune system itself. Together, these data point to an autoimmunity-prone state in DS, in which a steady-state cytokinopathy, hyperactivated CD4 T cells and ongoing B cell activation all contribute to a breach in immune tolerance. Our findings also open therapeutic paths, as we demonstrate that T cell activation is resolved not only with broad immunosuppressants such as Jak inhibitors, but also with the more tailored approach of IL-6 inhibition.

Unambiguous detection of SARS-CoV-2 subgenomic mRNAs with single cell RNA sequencing.

Single cell RNA sequencing (scRNA-Seq) studies have provided critical insight into the pathogenesis of Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2), the causative agent of COronaVIrus Disease 2019 (COVID-19). scRNA-Seq workflows are generally designed for the detection and quantification of eukaryotic host mRNAs and not viral RNAs. Here, we compare different scRNA-Seq methods for their ability to quantify and detect SARS-CoV-2 RNAs with a focus on subgenomic mRNAs (sgmRNAs). We present a data processing strategy, single cell CoronaVirus sequencing (scCoVseq), which quantifies reads unambiguously assigned to sgmRNAs or genomic RNA (gRNA). Compared to standard 10X Genomics Chromium Next GEM Single Cell 3' (10X 3') and Chromium Next GEM Single Cell V(D)J (10X 5') sequencing, we find that 10X 5' with an extended read 1 (R1) sequencing strategy maximizes the detection of sgmRNAs by increasing the number of unambiguous reads spanning leader-sgmRNA junction sites. Using this method, we show that viral gene expression is highly correlated across cells suggesting a relatively consistent proportion of viral sgmRNA production throughout infection. Our method allows for quantification of coronavirus sgmRNA expression at single-cell resolution, and thereby supports high resolution studies of the dynamics of coronavirus RNA synthesis.

Inducible CRISPR activation screen for interferon-stimulated genes identifies OAS1 as a SARS-CoV-2 restriction factor.

Interferons establish an antiviral state through the induction of hundreds of interferon-stimulated genes (ISGs). The mechanisms and viral specificities for most ISGs remain incompletely understood. To enable high-throughput interrogation of ISG antiviral functions in pooled genetic screens while mitigating potentially confounding effects of endogenous interferon and antiproliferative/proapoptotic ISG activities, we adapted a CRISPR-activation (CRISPRa) system for inducible ISG expression in isogenic cell lines with and without the capacity to respond to interferons. We used this platform to screen for ISGs that restrict SARS-CoV-2. Results included ISGs previously described to restrict SARS-CoV-2 and novel candidate antiviral factors. We validated a subset of these by complementary CRISPRa and cDNA expression experiments. OAS1, a top-ranked hit across multiple screens, exhibited strong antiviral effects against SARS-CoV-2, which required OAS1 catalytic activity. These studies demonstrate a high-throughput approach to assess antiviral functions within the ISG repertoire, exemplified by identification of multiple SARS-CoV-2 restriction factors.

ISG15 deficiency restricts HIV-1 infection.

Type I interferons (IFN-Is) are a group of potent inflammatory and antiviral cytokines. They induce IFN stimulated genes (ISGs), which act as proinflammatory mediators, antiviral effectors, and negative regulators of the IFN-I signaling cascade itself. One such regulator is interferon stimulated gene 15 (ISG15). Humans with complete ISG15 deficiency express persistently elevated levels of ISGs, and consequently, exhibit broad spectrum resistance to viral infection. Here, we demonstrate that IFN-I primed fibroblasts derived from ISG15-deficient individuals are more resistant to infection with single-cycle HIV-1 compared to healthy control fibroblasts. Complementation with both wild-type (WT) ISG15 and ISG15ΔGG (incapable of ISGylation while retaining negative regulation activity) was sufficient to reverse this phenotype, restoring susceptibility to infection to levels comparable to WT cells. Furthermore, CRISPR-edited ISG15ko primary CD4+ T cells were less susceptible to HIV-1 infection compared to cells treated with non-targeting controls. Transcriptome analysis of these CRISPR-edited ISG15ko primary CD4+ T cells recapitulated the ISG signatures of ISG15 deficient patients. Taken together, we document that the increased broad-spectrum viral resistance in ISG15-deficiency also extends to HIV-1 and is driven by a combination of T-cell-specific ISGs, with both known and unknown functions, predicted to target HIV-1 replication at multiple steps.

The NF-κB Transcriptional Footprint Is Essential for SARS-CoV-2 Replication.

SARS-CoV-2, the etiological agent of COVID-19, is characterized by a delay in type I interferon (IFN-I)-mediated antiviral defenses alongside robust cytokine production. Here, we investigate the underlying molecular basis for this imbalance and implicate virus-mediated activation of NF-κB in the absence of other canonical IFN-I-related transcription factors. Epigenetic and single-cell transcriptomic analyses show a selective NF-κB signature that was most prominent in infected cells. Disruption of NF-κB signaling through the silencing of the NF-κB transcription factor p65 or p50 resulted in loss of virus replication that was rescued upon reconstitution. These findings could be further corroborated with the use of NF-κB inhibitors, which reduced SARS-CoV-2 replication in vitro. These data suggest that the robust cytokine production in response to SARS-CoV-2, despite a diminished IFN-I response, is the product of a dependency on NF-κB for viral replication. IMPORTANCE The COVID-19 pandemic has caused significant mortality and morbidity around the world. Although effective vaccines have been developed, large parts of the world remain unvaccinated while new SARS-CoV-2 variants keep emerging. Furthermore, despite extensive efforts and large-scale drug screenings, no fully effective antiviral treatment options have been discovered yet. Therefore, it is of the utmost importance to gain a better understanding of essential factors driving SARS-CoV-2 replication to be able to develop novel approaches to target SARS-CoV-2 biology.

Human TBK1 deficiency leads to autoinflammation driven by TNF-induced cell death.

TANK binding kinase 1 (TBK1) regulates IFN-I, NF-κB, and TNF-induced RIPK1-dependent cell death (RCD). In mice, biallelic loss of TBK1 is embryonically lethal. We discovered four humans, ages 32, 26, 7, and 8 from three unrelated consanguineous families with homozygous loss-of-function mutations in TBK1. All four patients suffer from chronic and systemic autoinflammation, but not severe viral infections. We demonstrate that TBK1 loss results in hypomorphic but sufficient IFN-I induction via RIG-I/MDA5, while the system retains near intact IL-6 induction through NF-κB. Autoinflammation is driven by TNF-induced RCD as patient-derived fibroblasts experienced higher rates of necroptosis in vitro, and CC3 was elevated in peripheral blood ex vivo. Treatment with anti-TNF dampened the baseline circulating inflammatory profile and ameliorated the clinical condition in vivo. These findings highlight the plasticity of the IFN-I response and underscore a cardinal role for TBK1 in the regulation of RCD.

Transcriptional profiling of mESC-derived tendon and fibrocartilage cell fate switch.

The transcriptional regulators underlying induction and differentiation of dense connective tissues such as tendon and related fibrocartilaginous tissues (meniscus and annulus fibrosus) remain largely unknown. Using an iterative approach informed by developmental cues and single cell RNA sequencing (scRNA-seq), we establish directed differentiation models to generate tendon and fibrocartilage cells from mouse embryonic stem cells (mESCs) by activation of TGFß and hedgehog pathways, achieving 90% induction efficiency. Transcriptional signatures of the mESC-derived cells recapitulate embryonic tendon and fibrocartilage signatures from the mouse tail. scRNA-seq further identify retinoic acid signaling as a critical regulator of cell fate switch between TGFß-induced tendon and fibrocartilage lineages. Trajectory analysis by RNA sequencing define transcriptional modules underlying tendon and fibrocartilage fate induction and identify molecules associated with lineage-specific differentiation. Finally, we successfully generate 3-dimensional engineered tissues using these differentiation protocols and show activation of mechanotransduction markers with dynamic tensile loading. These findings provide a serum-free approach to generate tendon and fibrocartilage cells and tissues at high efficiency for modeling development and disease.

An ACE2 Triple Decoy that neutralizes SARS-CoV-2 shows enhanced affinity for virus variants.

The SARS-CoV-2 variants replacing the first wave strain pose an increased threat by their potential ability to escape pre-existing humoral protection. An angiotensin converting enzyme 2 (ACE2) decoy that competes with endogenous ACE2 for binding of the SARS-CoV-2 spike receptor binding domain (S RBD) and inhibits infection may offer a therapeutic option with sustained efficacy against variants. Here, we used Molecular Dynamics (MD) simulation to predict ACE2 sequence substitutions that might increase its affinity for S RBD and screened candidate ACE2 decoys in vitro. The lead ACE2(T27Y/H34A)-IgG1FC fusion protein with enhanced S RBD affinity shows greater live SARS-CoV-2 virus neutralization capability than wild type ACE2. MD simulation was used to predict the effects of S RBD variant mutations on decoy affinity that was then confirmed by testing of an ACE2 Triple Decoy that included an additional enzyme activity-deactivating H374N substitution against mutated S RBD. The ACE2 Triple Decoy maintains high affinity for mutated S RBD, displays enhanced affinity for S RBD N501Y or L452R, and has the highest affinity for S RBD with both E484K and N501Y mutations, making it a viable therapeutic option for the prevention or treatment of SARS-CoV-2 infection with a high likelihood of efficacy against variants.

Atypical Inflammatory Syndrome Triggered by SARS-CoV-2 in Infants with Down Syndrome.

While adults with Down syndrome (DS) are at increased risk of severe COVID-19 pneumonia, little is known about COVID-19 in children with DS. In children without DS, SARS-CoV-2 can rarely cause severe COVID-19 pneumonia, or an even rarer and more typically pediatric condition, multisystem inflammatory syndrome in children (MIS-C). Although the underlying mechanisms are still unknown, MIS-C is thought to be primarily immune-mediated. Here, we describe an atypical, severe form of MIS-C in two infant girls with DS who were hospitalized for over 4 months. Immunological evaluation revealed pronounced neutrophilia, B cell depletion, increased circulating IL-6 and IL-8, and elevated markers of immune activation ICAM1 and FcÉ£RI. Importantly, uninfected children with DS presented with similar but less stark immune features at steady state, possibly explaining risk of further uncontrolled inflammation following SARS-CoV-2 infection. Overall, a severe, atypical form of MIS-C may occur in children with DS.

Pathogenesis, MicroRNA-122 Gene-Regulation, and Protective Immune Responses After Acute Equine Hepacivirus Infection.

BACKGROUND AND AIMS: Equine hepacivirus (EqHV) is phylogenetically the closest relative of HCV and shares genome organization, hepatotropism, transient or persistent infection outcome, and the ability to cause hepatitis. Thus, EqHV studies are important to understand equine liver disease and further as an outbred surrogate animal model for HCV pathogenesis and protective immune responses. Here, we aimed to characterize the course of EqHV infection and associated protective immune responses. APPROACH AND RESULTS: Seven horses were experimentally inoculated with EqHV, monitored for 6 months, and rechallenged with the same and, subsequently, a heterologous EqHV. Clearance was the primary outcome (6 of 7) and was associated with subclinical hepatitis characterized by lymphocytic infiltrate and individual hepatocyte necrosis. Seroconversion was delayed and antibody titers waned slowly. Clearance of primary infection conferred nonsterilizing immunity, resulting in shortened duration of viremia after rechallenge. Peripheral blood mononuclear cell responses in horses were minimal, although EqHV-specific T cells were identified. Additionally, an interferon-stimulated gene signature was detected in the liver during EqHV infection, similar to acute HCV in humans. EqHV, as HCV, is stimulated by direct binding of the liver-specific microRNA (miR), miR-122. Interestingly, we found that EqHV infection sequesters enough miR-122 to functionally affect gene regulation in the liver. This RNA-based mechanism thus could have consequences for pathology. CONCLUSIONS: EqHV infection in horses typically has an acute resolving course, and the protective immune response lasts for at least a year and broadly attenuates subsequent infections. This could have important implications to achieve the primary goal of an HCV vaccine; to prevent chronicity while accepting acute resolving infection after virus exposure.

Single-cell resolution landscape of equine peripheral blood mononuclear cells reveals diverse cell types including T-bet+ B cells.

BACKGROUND: Traditional laboratory model organisms represent a small fraction of the diversity of multicellular life, and findings in any given experimental model often do not translate to other species. Immunology research in non-traditional model organisms can be advantageous or even necessary, such as when studying host-pathogen interactions. However, such research presents multiple challenges, many stemming from an incomplete understanding of potentially species-specific immune cell types, frequencies, and phenotypes. Identifying and characterizing immune cells in such organisms is frequently limited by the availability of species-reactive immunophenotyping reagents for flow cytometry, and insufficient prior knowledge of cell type-defining markers. RESULTS: Here, we demonstrate the utility of single-cell RNA sequencing (scRNA-Seq) to characterize immune cells for which traditional experimental tools are limited. Specifically, we used scRNA-Seq to comprehensively define the cellular diversity of equine peripheral blood mononuclear cells (PBMC) from healthy horses across different breeds, ages, and sexes. We identified 30 cell type clusters partitioned into five major populations: monocytes/dendritic cells, B cells, CD3+PRF1+ lymphocytes, CD3+PRF1- lymphocytes, and basophils. Comparative analyses revealed many cell populations analogous to human PBMC, including transcriptionally heterogeneous monocytes and distinct dendritic cell subsets (cDC1, cDC2, plasmacytoid DC). Remarkably, we found that a majority of the equine peripheral B cell compartment is comprised of T-bet+ B cells, an immune cell subpopulation typically associated with chronic infection and inflammation in human and mouse. CONCLUSIONS: Taken together, our results demonstrate the potential of scRNA-Seq for cellular analyses in non-traditional model organisms and form the basis for an immune cell atlas of horse peripheral blood.

ISG15 protects human Tregs from interferon alpha-induced contraction in a cell-intrinsic fashion.

OBJECTIVES: Type I interferons (IFNs) inhibit regulatory T-cell (Treg) expansion and activation, making them beneficial in antiviral responses, but detrimental in autoimmune diseases. Herein, we investigate the role of ISG15 in human Tregs in the context of refractoriness to type I IFN stimulation. METHODS: ISG15 expression and Treg dynamics were analysed in vitro and ex vivo from patients with chronic hepatitis C, with lupus and ISG15 deficiency. RESULTS: ISG15 is expressed at high levels in human Tregs, renders them refractory to the IFN-STAT1 signal, and protects them from IFN-driven contraction. In vitro, Tregs from healthy controls upregulate ISG15 upon activation to higher levels than conventional CD4 T cells, and ISG15-silenced Tregs are more susceptible to IFNα-induced contraction. In human ISG15 deficiency, patient Tregs display an elevated IFN signature relative to Tregs from healthy control. In vivo, in patients with chronic hepatitis C, 2 days after starting pegIFN/ribavirin therapy, a stronger ISG15 inducibility correlates with a milder Treg depletion. Ex vivo, in systemic lupus erythematosus patients, higher levels of ISG15 are associated to reduced STAT1 phosphorylation in response to IFNα, and also to increased frequencies of Tregs, characterising active disease. CONCLUSION: Our results reveal a Treg-intrinsic role of ISG15 in dictating their refractoriness to the IFN signal, thus preserving the Treg population under inflammatory conditions.

Single-cell RNA sequencing of equine mesenchymal stromal cells from primary donor-matched tissue sources reveals functional heterogeneity in immune modulation and cell motility.

BACKGROUND: The efficacy of mesenchymal stromal cell (MSC) therapy is thought to depend on the intrinsic heterogeneity of MSC cultures isolated from different tissue sources as well as individual MSCs isolated from the same tissue source, neither of which is well understood. To study this, we used MSC cultures isolated from horses. The horse is recognized as a physiologically relevant large animal model appropriate for translational MSC studies. Moreover, due to its large size the horse allows for the simultaneous collection of adequate samples from multiple tissues of the same animal, and thus, for the unique collection of donor matched MSC cultures from different sources. The latter is much more challenging in mice and humans due to body size and ethical constraints, respectively. METHODS: In the present study, we performed single-cell RNA sequencing (scRNA-seq) on primary equine MSCs that were collected from three donor-matched tissue sources; adipose tissue (AT), bone marrow (BM), and peripheral blood (PB). Based on transcriptional differences detected with scRNA-seq, we performed functional experiments to examine motility and immune regulatory function in distinct MSC populations. RESULTS: We observed both inter- and intra-source heterogeneity across the three sources of equine MSCs. Functional experiments demonstrated that transcriptional differences correspond with phenotypic variance in cellular motility and immune regulatory function. Specifically, we found that (i) differential expression of junctional adhesion molecule 2 (JAM2) between MSC cultures from the three donor-matched tissue sources translated into altered cell motility of BM-derived MSCs when RNA interference was used to knock down this gene, and (ii) differences in C-X-C motif chemokine ligand 6 (CXCL6) expression in clonal MSC lines derived from the same tissue source correlated with the chemoattractive capacity of PB-derived MSCs. CONCLUSIONS: Ultimately, these findings will enhance our understanding of MSC heterogeneity and will lead to improvements in the therapeutic potential of MSCs, accelerating the transition from bench to bedside.